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Abstract 

Following the statistical mechanics methodology, firstly introduced in macroeco- 
nomics by Aoki [1, 2, 3], we provide some insights to the well known works of 
(author?) [7, 6]. Specifically, we reach analytically a closed form solution of their 
models overcoming the aggregation problem. The key idea is to represent the econ- 
omy as an evolving complex system, composed by heterogeneous interacting agents, 
that can partitioned into a space of macroscopic states. This meso level of ag- 
gregation permits to adopt mean field interaction modeling and master equation 
techniques. 
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1 Introduction 



Aim of the present work is to replicate the mechanism of the models intro- 
duced in j^j and [i| and to resettle them in a dynamic stochastic framework, 
as defined by Aoki [l], 0, Hf . The quantitative previsions of j3] (and of its de- 
velopments [12j]) replicate a number of stylized facts, strengthening the idea 
that the economy would be better represented as a complex dynamical sys- 
tem rather than a mere sum of identical and perfect informed agents. Any 
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attempt to resolve analytically this kind of models must face the aggregation 
issue, since the problem of how to sum up heterogeneous and evolving agents 



cannot be dealt with the consuete tools of the economist 14]. Masanao Aoki 
[USS] managed to remove the representative agent hypothesis, introducing in 
economics the ground-breaking concept of mean-field interaction, that makes 
feasible the analytical aggregation of heterogeneous agents, replacing the un- 
realistic mechanic determinism of mainstream framework with a set stochastic 
tools borrowed from statistical mechanics. Mean-field interaction can be de- 
fined as the average interaction model that substitutes all the relations among 
agents that, otherwise, could not be analytically treated[5|. All the agents 
are clusterized in a pre-defined set of states, basing on one particular feature 
(the micro-state, e.g. the level of production for each firm) that determines 
the characteristics and the evolution of the aggregate (the macro-state, e.g. 
the total level of output). The accent is not on the single agent, but on the 
number or proportion of agents that occupy a certain state at a certain time. 
These levels are governed by a stochastic law, that defines also the functional 
of the probability distributions of aggregate variables and, if existing, their 
equilibrium distributions. The structure of the work is the following: firstly, 
we specify the hypothesis for the stochastic structure of the system (section 
|5]and for the firms that compose it (section [3])); then (section HJ), we develop 
the probabilistic framework, setting the dynamical instruments needed for the 
aggregation. 



2 Structure of the system and definition of states 



We set up a model in continuous time for a system of heterogeneous and inter- 
acting agents, partitioned into groups or states. In this paragraph we state the 
hypothesis that are at the root of this stochastic dynamics, and in particular 
the definition and the structure of the states. Our system is articulated in two 
states. This partition permits to isolate the effect of bankruptcy costs on the 
aggregate dynamics. Moreover, the application of recently proposed methods 
for systems with higher order of states [20(] would entail a too high price in 
terms of computational complication respect to the expected improvement in 
the realism of the model. Along time, a single firm can be in one of the two 
states, depending on its financial soundness, proxied by the equity ratio (the 
ratio among the net worth and the total assets). Therefore there are two types 
of firms: the "good" firms, that have a high equity ratio, and the "bad" firms, 
that have a low equity ratio that exposes them to the risk of demisd"*"! System 



works in count inuos time [10|: t G T C K+. 

The economy is populated by a fixed number of firms N = N(t) : Vt G T, 
each indexed by % for any given time. 



1 To be complete, there is also a third type, the "ugly" firms, i.e. the failed ones, 
as explained in the following. 
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The system's vector of states lj is identified by the financial condition of the 
firm: 



u{t) = {ui{t) = HV(ai(t)\ a) Vz <N} : HV( ai (t)\a) = < 



1 cn(t) < a 

ai(t) > a 



where a represents the threshold of equity ratio, that individuates firms that 
are in a critical financial situation, and, therefore, for which the probability of 
bankruptcy is bigger than and where aiit) is the equity ratio for firm i at 
time t. 

In what follows we set ai(t) = x for firms which equity ratio is under the 
threshold and ci2 (t) = y for firms which equity ratio is over the threshold. 
The cardinality of the j-th state, to say the number of firms in state j — 0,1 
is given by 



card N^t) = #{uji(t) = 1 V i e 1} = N^t) 
N (t) — N — JVi (t) 



N(t) = (N (t),N 1 (t))(l) 



By assumption, the dynamics of the occupation number Nj follows a continu- 
ous time jump Markov process, defined over a state space Q = (x, y) equipped 
with the following counting measure Nq(.) : Q x T — > N, so that (£7,iV()(.)) 
is a countable sample space. In the following we indicate with x the case of 
firms with equity ratio below the threshold a and with y the alternative case: 

uj = x <^> Ui(t) = 1 V uo = y <^> u)i(t) = (2) 



such that Nf\(.) evaluates the cardinality of microstates: N^(t) = Ni(t) 
uj = x and JV( W ) = N (t) ^=^> u = y. The relative frequency of firms 
is indicated in small letters: rik = N^/N. A-priori probability of u> — 1 is 
indicated by r\: 

p(cu = x) = r] <^ p{uj = y) = 1 — rj 



A firm entries into the system in state x and fails (and thus exits from the 
system) only if it is in state x (i. e. if its probability of bankruptcy is bigger 
than 0); in order to maintain constant the number of firms iV we assume that 
each bankrupted firm is immediately substituted by a new one. Therefore 
failures of firms do not modify the value of JVi := N(x). Firms move from x 
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to y or vice versa according to the following transition rates: 



b(N 3 ) = r(N 3 + l\N 3 ) = (^(1 - v) 
d(Nj) = l(Nj - l\Nj) = l^, 



(3) 



where ( is the transition probability from state y to x (firms whose finan- 
cial position is deteriorated from a period to another, with equity ratio that 
becomes lower than a) and i the probability of the inverse transition (firms 
whose equity ratio improved becoming bigger than a). Having already indi- 
cated with iV- 5 the occupation number of firms in state j, transition rates can 
be evaluated according to: 

b(Nj) = r(N 3 + l\N 3 ) = X(N - N 3 ) : X — ((l — rj) 
d(N 3 )=l(N j -l\N 3 )= 1 (N 3 ): 1 = i V 



In statistical mechanics terms, this kind of system can be defined as a statis- 
tical ensemble with conservative cardinality, described by a continuous time 
Markov process over a discrete state space with the structure of a birth-death 
process. 



3 The Firms 



3. 1 Hypothesis 



The assumptions regarding firms are the same of the original models of I (author?) 
(cjl, 0], as modified by 0|, except for the following specifications. Respect to 
the theoretical construction proposed by [3] this model does not consider any 
credit market. There is no direct but a mean-field interaction. Differently from 
(H the total number of firms N does not change over time. Firms are identical 
within each state. The production function of a generic firm i is: 

<s(t) = 2(k t (t)) 1 / 2 (5) 



where k is the physical capital and q is the physical output. It follows that the 
demand of capital function for a single firms will be equal to: 

h(t) = (6) 
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The multiplicative shock price Ui(t) has uniform distribution with support 
[0.75; 1.25] and E(ft) = 1. Given all the above specified hypothesis, the profit 
function for a generic firm % can be expressed by: 

TTi(t) = P{t)ui{t) qi {t) - ri(t)Ki(t) (7) 



Once a firm get failed, it faces bankruptcy costs Cj(t)/i growing with the size 
of firm, proxied by the volume of its sales, and quantified by: 

d(t) = ci^P^qit),) 2 = c(P(t) Ul (t) qi (t)) 2 : 0<c<l (8) 



3.1.1 Transition probabilities 

A firm goes bankrupted if it "consumes" all its own capital, then < 0. 
So that, analogously with [ij], we can express the bankruptcy condition as a 
function of the price shock: 

a(0<(||)(rJf<W/«W-*W^| 5 ) S ft(0 



Substituting equation [6] into the above expression and, without loss of gener- 
ality, normalizing reference price P(t) = P(t') to 1, it is possible to simplify 
the r.h.s. of the above equation: 

u(t') = ^(r-a x (t)) (9) 



Recalling that the random variable u has support [0.75; 1.25], the critical 
thresholds of shock prices for having bankruptcy will be: 

u = 0.75 <^=> Ui(t) < 0.75 
< u G (0.75; 1.25) 0.75 < u^t) < 1.25 (10) 

u = 1.25 Ui(t) > 1.25 



Then, it is possible to indicate the probability fi(t) of failure for a firm as a 
function of u{t) : 

p(t) = F(u(t)) = p(u(t) < u(t)) = M(t) ~°- 75 = 2u(t) - 1.5 (11) 

0.5 



Equation (fl2"l) permits to determine endogenously the threshold a: indeed it 
can be interpreted as the minimum value of equity ratio which ensures the 
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surviving of the firm, i. e. for which the probability of bankruptcy is equal to 
zero, and, therefore, it can be expressed as: 

^- r -W) (12) 



With an analogue procedure, we can specify the transition probabilities ( and 
i as dependent variables of the price shock Ui(t). Indicating the critical values, 
respectively with Ut^(t) and u i)lr (t) we obtain: 



Ui(t)<3f-(r + a(t)-a (t))=u c (t) 
u l (t)> 3l f(r + a(t)-a 1 (t))=%(t) 



since the thresholds, that for the [TT] is equal to 0, here become equal to, 
respectively (a(t) — a (t)) and (cti(t) — a(t)). The range of variation of the 
two thresholds is truncated as done in [TDJ It is straightforward now to get the 
transition probability for each state: 

C(t) = p(u{t) < u c {t)) = 2u ( {t) - 1.5 (13) 
L (t) = 1 - p{u(t) < u b (t)) = -2u L (t) + 2.5 (14) 

3.2 Firms object function. 



A firm decide the optimal quantity to produce in order to maximize its profit, 
using the information at its disposal. Under the stated hypothesis the object 
function of a generic firm i can be then expressed as: 

maxF(g,(t)) := {E [P(t) Ui (t') qi (t)} - rK^t) - C^t)^')} (15) 

We suppose that firms take into consideration the present level of failure prob- 
ability, therefore that /i(t') = /i(t), and that P e {t') = P{t') = P(t) = 1 without 
loss of generality. The first order condition is: 

1 - r qi (t) - 2cqi(t)fi(t) = 



Consequently we get two different optimal levels of production for firms in 
state x and for firms in state y, respectively: 

q{ = (r + 2c/i)- 1 (16) 

* -l 
q 2 =r 
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since ji = for firms in state y. 



4 Stochastic inference 



4-1 Master equation and stationary points 



In order to specify the dynamics of the joint probabilities, and, by this way, 
the stochasticH] evolution of the system, following 0, chap. 3], [8, page 252] 



and 11, page 62] we make use of the following master equation: 



dP(N k ,t) 



dt 



6(AA fe _ 1 )P(iV i: _ 1 ) + d(N k+1 )P(N k+1 )+ 
{[(b(N k ) + d(N k ))P(N k )}} 



(17) 



with boundary conditions: 

P(N, t) = 6(iV 1 )P(A^ 1 - 1, t) + d(N)P(N, t) 
P(0,t) = 6(l)P(l,t) + d(0)P(0,t) 



These conditions ensure that the distributions functions consider only con- 
sistent values, that is to say Aq G [0; N] and, therefore, n\ G [0;1]. Since 
an analytical solution for master equations can be obtained only under very 



specific conditions [231] , we solve it with an approximation method, based on 
led and lag operatorajj. We assume that the number of firms in state x at a 
given moment is determined by its expected mean (m), the drift, and by an 
additive fluctuations component of order N~ l l 2 around this value, that is to 
say the spread: 

N l = Nm + VNs (19) 



Solving the modified master equation (see Appendix a) we find an explicit 
formulation for the macroscopic equation: 

^ = -(A + 7 )iV 1 + AiV (20) 



2 Indeed, "the master equation describes the time evolution of the probability dis- 
tribution of states, not the time evolution of the states themselves" 0, page 7]. 

3 For a detailed exposition see 0] or 
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This equation describes the dynamics of the drift. It can be interpretated as 
the long-run trend of the occupation number and, then, keeping all the other 
relevant variables unchanged, of the production. Now we can determine the 
stationary value of E(n fc ) = m and the consequent steady state equilibrium of 
the economy, simply setting the r.h.s. of the (1201) to 0: 

A + 7 



Once demonstrated the existence of a dynamical equilibrium point, a deeper 
insight on system's dynamic is needed in order to verify the convergence of the 
dynamics toward the equilibrium distribution. The solution of the differential 
equation (1201) is: 



m (t) = m e- (A+7) ' + — ^— (22) 

A + 7 



that, setting an initial point m , yields to: 

m(t) = m* + (m - m*)e~ (A+7) * (23) 



that verifies the convergence and the stability of the equilibrium since the 
second term goes to as t — > oo. 



4-2 Aggregate output 



The aggregate output of the system will be: 

K W = _* (24) 
r + ZCjJL(t) r 



Keeping constant all the other variables, dynamic fluctuations in the level of 
Y(t) are due to the changing in the levels of Ni and N . Once quantified the 
equilibrium distribution of the drift, we are able now to obtain the steady 
state value of aggregate production: 



Y e = N 



A 



2c/i 



r A + 7 r(r + 2c/i) 



N 



A 



A + 7 



(25) 



As can be easily seen the factor inside the brackets of the (1251) is equal to 
[y - m*(y 1 - y )] 
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that testifies the consistency of the model. Since the aggregate production 
function depends on m, also its dynamic will be convergent to a stationary 
level. Dynamics comes out to be dependent on the transition rates A and 7 and 
on the differences in firms level of production. The dynamics of these factors 
are studied in the next section. 



4-3 Equilibrium distribution and critical points 



Equating the master equation to 0, it is possible to obtain the Kolmogorov 
condition that equates the probability fluxes entering in a state with the fluxes 
coming out from that state. Formally, using equations H] and [20) 

m = Cv( m ) ~ [^(1 ~~ r](m)) + Cv( m )] m (26) 



Setting m to and rearranging, the stationary configuration of the system can 
be expressed as: 

77(771*) Lm* , , 

m = => , . = — (27) 

1 — r](m*) (,(1 — m*) 



Being in Markovian space, we can make use of Brook's lemma [l5| that defines 
local characteristic for this kind of chains: 



By means of Hammersley and Clifford theorem (see [16J as demonstrated in 
[171]). the stationary probability of the markovian process for Nj, when detailed 
balance holds, can be expressed by: 

P e {N,) oc Z~ l e- pNU ^ (29) 



where U(N^) is the Gibbs potential [9j. The parameter (3 may be interpreted 
as an inverse measure of the uncertainty of the system. The above formulation 
allows us to express explicitly the values of the a-priori probabilities: 



1 - n{N 3 ) 

so that: 



rj(Nj) = N- X e p9 ^ 



+ e 



N 



(30) 
(31) 



(32) 
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It is easy to verify that large values of (3 associated with positive values of 
g{Nj) cause t](Nj) to be larger than 1 — i](Nj), making transition from state 
y to state x more likely to occur than the opposite one. On the other hand, 
values of (3 close to 0, make i](Nj) close to 0.5. To get a deeper insight on the 
meanings of (3 let us carry out the potential equation. In binary models and 
for great N, the equation of the potential is: 

A" 

U(N*) = -2 J g{N 3 )dy - -H(N) 
o P 



where H(N) is the Shannon entropy for the vector of occupation numbers. 
g{Nj) is a function that evaluates the relative difference in the outcome as a 
function of Nj. In intuitive terms, to individuate the stationary points of prob- 
ability dynamics we need to individuate its peak. (3 is an inverse multiplicative 
factor for entropy: this implies then that, for very large values of /3, the entropy 
component does not play any role. On the contrary, as (3 approaches 0, the 
weight of the entropy component grows. In terms of the present model, a rela- 
tive high value of (3 means that the uncertainty in the system is low, with few 
firms exposed at bankruptcy risk and, due to this, the bankruptcy probability 
\x plays a negligible role pp. 55 and following]. Therefore for values of f3 
around 0, and a more relevant volatility in the system, in order to individuate 
the peak of probability dynamic we need to find the local minimum of the 
potential. [2| shows that the points in which the potential is minimized are 
also the critical point of the aggregate dynamics of pj. Deriving the potential 
respect to iV J : 

q(N i) - _±^mi - _ J_ ln I m 



Solving the following MaxEnt problem [21fl : 



max#(JVi,iVo) = -N^t) ln(iVi(t)) - N (t) \n(N (t)) s.t. 
N 1 (t)+N (t) = N (34) 
N 1 (t)y 1 {t) + N (t)y (t) = Y(t) 



with suitable Lagrange multipliers equal to, respectively, Ai = 1 — a and 
A 2 = — (3 we get a closed solution for (3: 

(3(t) = In (- y -^M ) ( yi (t) - yoit))- 1 (35) 

V yo{t)-y[t)J 
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Setting U' = and using equation (I35I) . we get an explicit formulation for 
g(Nx) in stationary condition: 



9(Ni) = VJL ^- (36) 



that quantifies the mean difference (for states) of the outcome. 

From equation [32] it follows that the point of local minimal of the potential is 

given by: 

U' = e 2/ W) = Nj (37) 
N — Nj v ' 



that, if the rates of entries and exits are equated, (i. e. if i = (), reproduces 
exactly the [271 Therefore, making use of the [30] we can write: 

e Wg{Nj) 



that is the maximum likelihood estimation of the Gibbs distribution of the 
firms in state Nj. Let us analyze the different behavior of the stationary dis- 
tribution for different values of (3. 

Given equation [33J the critical points in which the potential is 



(3 — > oo 



minimized are also the zeros of the function g(Nj): 

U'(m*) = -2g{m*) = (39) 

This suggests that (3 may be interpreted as an inverse index of uncertainty [it]. 
Considering equation | 



g (Nj) =0^2/o-2/i = 

Under the specified conditions, there is no uncertainty in the system, since no 
firm can go bankrupted and, therefore, the level of production is unambigu- 
ously determined. Indeed the value of (3 can go to infinity if iVo — > iV or if 
// — > 0, since both situation imply a convergence among the different targets 
of production at micro level and, then, a minimum degree of uncertainty in 
the system. 

(3 — > In this second case, in order to individuate the critical points of 



the dynamics, a further deepening is needed. (3 can go to if and only if 
J^ N . — ¥ 1, that is to say, if the system is populated, quite in the same pro- 
portion, of firms in state x and of firms in state y. But this is not informative 
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about the behavior of the g(N) since the first factor of the equation [331 goes 
to infinity (given that (3 — > 0) while the second goes to 0. The Cox and Miller 
hazard function [l8|, with a density function analogous to the one of equal to a 
standard Brownian motion's first-passage [191]. in terms of m can be expressed 



as: 



F(m) 



1 + e 



-2/3m 



h(m) 



2(3 



1 + e- 2 ^ m 



(40) 



Now, we have to calculate the probability that a firm passes from a state to 
another in response to a small variation in the difference of relative production, 
conditional on the current difference among y and y\ (quantified by g(m), 
from equation I3"oT) . We can then rewrite the conditional hazard function in the 
following way: 



P{v <y%- y \m* 



l + e 



-Wg{m") 



(41) 



and then: 



Kyi - yo\m*) 



2(3r](m* 



1 _|_ e -2/3g(m*) 



(42) 



Supposing that 77 (m*) = m*, we finally obtain: 
h{yi - Vol™*) = 2pm* 



(43) 



Therefore, we may conclude that for values of (3 close to 0, the critical point 
of probability dynamics, here found by minimizing the potential, is a value of 
m* approximately equal to (3 itself. In other words, (3 may be considered the 
as the conditional hazard rate in the range where (3 is small. The potential 
then is minimized for a fraction m* of firms in state x when the value of the 
conditional hazard function is approximately equal to (3. 



5 Concluding remarks 



In this work we apply Aoki's methodology to the model presented in [6[ and 
modified, introducing interacting heterogeneous agents, in Q]. The adoption 
of the statistical mechanical approach, based on mean field interaction and 
aggregate stochastic dynamic analysis, leads to a stable analytical solution, 
without recurring to computer simulations. We start from very general con- 
ditions and hypothesis, embodied in a stochastic framework, and, finally, we 
obtain an explicit macroscopic equation that describes the evolution of the 
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system and its long term equilibrium solution. The dynamics is analyzed by 
means of master equation solution techniques enriched by the use of MaxEnt 
and hazard function analysis. 



Appendix A 

The master equation has now to be expressed as Q(s), a function of s, and 
then becomes: 

iff at \ dQ dsdQ • , . 



with transition rates reformulated in the following way: 



b(s) = A 
d{s 



N - Nm 



Ns 



7 



Nm + VNs 



Since 



ds 
~dt 



dm 
~dt 



:-2) 

(-3) 



:-4) 



equation (J7TJ) can be expressed as: 

Q { s) = — -N'—m (.5) 



We rewrite again now the master equation ( TTfl) and the transition rates by 
means of lead and lag operators. The use of these operators allows to express 
the master equation in a more treatable form, making the two probability 
flows (in and out) homogeneous. Specifically the transition probabilities PJ 
become: 



L[d(N k )P(N k , t)) = d(N k+1 )P(N k+1 ) (.6) 

L- 1 [b(N k )P(N k , t)} = d(N k ^)P(N k ^) (.7) 

so that the master equation will be expressed in this way: 

Q(s) = (L — l)[d(s)Q(s)} + (L- 1 - l)[d(s)Q(s)} (.8) 



Using the modified transition rates (Tj6]) and expanding the so obtained master 
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equation in inverse powers of s to the second order we get: 

at-I dQ _ j^-l/2dmdQ_ _ 
dr dr ds 

N-W (£) [d(s)Q(s)} + iV-U (£) a [d(s)Q(s)}+ 

-N-l/2 (£) [b(s)Q(s)} + iV-i| (f ) 2 [b(s)Q(s)] + ... 

= N- 1 '* (&) [(d(s) - b(s))Q(s)] + iV-U ( * ) [(6(5) + + ... 



At this point, in order to match the component of the same orders of powers 
of N between and equations (p]) and (p9j), we need to rescale the variable 
t = tN. Knowing that: 



d(s) - b(s) = (A + 7 )(iVm + VNs) - \N = (X + ^)N k - XN 
d( s ) + b(s) = (A - 7 )(iVm + t/Ns) + XN = (X- ^)N k + XN 



and taking the derivatives, up to the second order, it is possible to obtain 
what Aoki defines as diffusion approximation [if: 

AT-l 9Q _ ]y-l/2dmdQ_ _ 
dr dr ds 

(A + i)Q( 8 ) + N-W{d( 8 ) - b(s)) (£) Q(s) + N- l l(b(s) + d(s)) (£) Q(s 

Equating the terms of order N~ l l 2 for the polynomial identity principle we 
get: 

= -^(M --(.)) (|)«(.) 
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